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Abstract: A spherical wrist of the serial type is said to be isotropic if it can attain a posture 
whereby the singular values of its Jacobian matrix are all identical and nonzero. What isotropy 
brings about is robustness to manufacturing, assembly, and measurement errors, thereby guaran- 
teeing a maximum orientation accuracy. In this paper we investigate the existence of redundant 
isotropic architectures, which should add to the dexterity of the wrist under design by virtue 
of its extra degree of freedom. The problem formulation leads to a system of eight quadratic 
equations with eight unknowns. The Bezout number of this system is thus 2 8 = 256, its BKK 
bound being 192. However, the actual number of solutions is shown to be 32. We list all solu- 
tions of the foregoing algebraic problem. All these solutions are real, but distinct solutions do 
not necessarily lead to distinct manipulators. Upon discarding those algebraic solutions that 
yield no new wrists, we end up with exactly eight distinct architectures, the eight corresponding 
(■■<■) ' manipulators being displayed at their isotropic posture. 

m ' 1 Introduction 

The kinematic design of redundant spherical wrists under isotropy conditions is the subject of 
this paper. A manipulator is call isotropic if its Jacobian matrix can attain isotropic values on 
certain postures pQ . A matrix, in turn, is called isotropic if its singular values are all identical 
and nonzero. Furthermore, the matrix condition number can be defined as the ratio of its 
greatest to its smallest singular values [2]. Thus, isotropic matrices have a minimum condition 
number of unity. The kinematic structure of industrial manipulators are frequently decoupled 
into a positioning and an orientation submanipulator. The latter is designed with revolute joints 
whose axes intersect. However, when these three joints are coplanar, the manipulator becomes 
singular. As a means to cope with singularities, redundant wrists have been suggested [3]. An 
extensive bibliography on the design of spherical wrists can be found in 0]. 

Prior to our analysis leading to the architectures sought, we recall a few geometric concepts 
in the section below. 
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2 Isotropic Sets of Points on the Unit Sphere 



Consider the set S = { Pk }" of n points on the unit sphere, of position vectors {efc}™. Ap- 
parently, all the vectors of the foregoing set are of unit Euclidean norm. The second-moment 
tensor H of S is defined as 

n 

H = £e fe e£ (1) 
1 

The set S is said to be isotropic if and only if its second- moment tensor is isotropic. Since H 
is symmetric and positive-definite, it is isotropic if its matrix representation is proportional to 
the 3x3 identity matrix 1, the proportionality factor, denoted here with a 2 , being the square 
of the triple singular value of H. In our case, apparently, the singular values of H coincide with 
its eigenvalues. 

We note that, if S is the set of vertices of a Platonic solid, then H is isotropic. Table [T] records 
the values of n and a for each Platonic solid. 





Tetrahedron 


Cube 


Octahedron 


Icosahedron 


Dodecahedron 


n 


4 


8 


6 


12 


20 


a 


4/3 


8/3 


2 


4 


20/3 



Table 1: The values of n and a for the Platonic solids 



Remark 1: It is apparent that, if a point Pk of an arbitrary set S of points on the unit sphere 
is replaced by its antipodal Qk, of position vector (\k = — e/j, then the second- moment tensor H 
of S is preserved. 

The replacement of a point on the unit sphere by its antipodal will be termed, henceforth, 
antipodal exchange. As a consequence of Remark 1, then, the isotropy of a set of points on the 
unit sphere is preserved under any antipodal exchanges. 

We started by recalling the second moment of a set of points on the unit sphere because this is 
simpler to handle than the corresponding first moment. Besides, in deriving isotropic spherical 
wrists, all we need is the second moment. The first moment of a set of points on the unit sphere 
is somewhat more elusive, because the centroid of the set must be a point on the unit sphere 
as well. Thus, while the second moment H was taken with respect to the center of the sphere, 
the first moment, when taken with respect to the centroid, must vanish. The centroid of the 
set not being of interest to us in the context of spherical-wrist design, it will be left out of the 
discussion. 

2.1 Trivial Isotropic Sets of Points on the Unit Sphere 

The simplest sets of isotropic points are thus the sets of vertices of the Platonic solids. Hence, 

Definition 1 (Trivial isotropic set) An isotropic setS ofn points on the unit sphere is called 
trivial if it consists of the set of vertices of a Platonic solid (inscribed, of course, on the unit 
sphere.) 

It is noteworthy that isotropic sets on the unit sphere exist that are none of the Platonic solids, 
e.g., the 64- vertex polyhedron defined by the molecule of the buckminsterfullercne, popularly 
known as the Buckyball. The name comes from the architect R. Buckminster Fuller, who used 
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this polyhedron as the structure of the geodesic dome built on occasion of the Universal Exhibit 
of 1967 in Montreal. This polyhedron is also present in the patterns of soccer balls. 

Also note that the set of points on the unit sphere leading to an isotropic architecture for 
a spherical wrist need not be laid out with the center of the sphere as its centroid, which is 
a condition found for points in the plane 5j. For example, the three points of intersection of 
the unit sphere with the axes of an orthogonal coordinate frame with origin at the center of 
the sphere define an isotropic spherical wrist, namely, the one most commonly encountered in 
commercial manipulators, yet the above set of points corresponds to none of the Platonic solids. 
This three-revolute wrist is termed orthogonal because its neighboring axes make right angles. 

Trivial sets of isotropic points are important because they allow the derivation of nontrivial 
sets by simple operations, as described below. 

2.2 Properties of Isotropic Sets of Points on the Unit Sphere 

First, note 

Lemma 1: An isotropic set S of points on the unit sphere remains isotropic under any isometric 
transformation of the set. 

An isometry being either a rigid-body rotation or a reflection, the foregoing lemma should be 
obvious. Moreover, rigid-body rotations of isotropic sets are uninteresting because they amount 
to looking at the given set from a different viewpoint. However, distinct isotropic sets can be 
derived from reflections of isotropic sets about planes or lines. Nevertheless, as shown in the 
Appendix, a reflection about a line amounts to a rigid-body rotation about the line through 
an angle of tt. As a consequence, then, only reflections about planes will be considered when 
defining new isotropic sets from trivial ones. 

2.3 Nontrivial Isotropic Sets of Points 

We show in this subsection, with a numerical example, how nontrivial sets of isotropic points 
on the unit sphere can be derived from a trivial set by application of reflections about planes. 
Now, since the reflection plane can be defined in infinitely-many ways, a correspondingly infinite 
number of reflections is possible. We are interested only in linearly-independent reflections, which 
are, apparently, only three, one about each of three mutually orthogonal planes. 

A trivially isotropic set of four points, namely, the vertices of a regular tetrahedron inscribed 
in the unit sphere, is given below: 





" 1 " 




" -1/3 " 




" -1/3 " 




' -1/3 " 


ei = 





j e 2 = 


-2V2/3 


, e 3 = 


V2/3 


, e 4 = 


y/2/3 














V6/3 




-V6/3 



We produce now three nontrivial sets of isotropic points by reflecting the foregoing set onto the 
three coordinate planes, y-z, x-z and x-y, successively. We display below the three new sets: 

1. The reflection with respect to the y-z plane gives 





" -1 " 




■ 1/3 " 




" 1/3 " 




■ 1/3 " 


ei = 





, e 2 = 


-2\/2/3 




V2/3 


, < = 


\/2/3 














V6/3 




-V6/3 



3 



2. The reflection with respect to the x-z plane gives 





' i " 




■ -1/3 - 




" -1/3 " 




" -1/3 " 







7 ^2 ~~ 


2V2/3 


. e' 3 ' = 


-a/2/3 


7 e 4 — 


-V2/3 














V6/3 




-76/3 



(2c) 



3. The reflection with respect to the x-y plane gives 





" 1 " 




■ -1/3 " 




r -1/3 " 




' -1/3 " 







p'" - 

7 ^2 — 


-2V2/3 


p'" - 

7 e 3 — 


V2/3 


p'" - 

7 e 4 — 


y/2/3 














-V6/3 




V6/3 



(2d) 



3 Isotropic Spherical Wrists 

Most serial wrists encountered in manipulators are provided with three revolute joints. We start 
with a general n-revolute spherical wrist, as depicted in Fig. [TJ with Jacobian matrix J given 
by© 



[ ei e 2 



(3) 



where, we recall, e& is the unit vector indicating the direction of 
the fcth revolute axis. We display below the kinematic relation 
between the joint-rate vector 6 and the angular-velocity vector 
u> of the end-effector (EE): 



36 = u> 



(4) 



It should be apparent that 

Rpmark 2: A set {e^}™ of unit vectors produces n! Jacobian 
matrices, and hence, n! distinct wrists. 

Kinetostatic isotropy requires that the singular values of the 
Jacobian matrix be all identical and nonzero, i.e., 




JJ T = cr 2 l 



Figure 1: A general n- 
(5) revolute spherical wrist 



where a is the common singular value, of multiplicity three, and 1 is, as defined earlier, the 3x3 
identity matrix. The isotropy condition thus leads to 



a 2 l 



(6) 



The value of a is found by taking the trace of both sides of eq.©, which yields 

n 

^2 e k ■ e/c = 3cr 2 



(7) 



and hence, 



i.e., if J is isotropic, then (a) every pair of n-dimensional rows of J is orthogonal and (b) the 
three rows of J have the same Euclidean norm, namely, \J n/3. 

Now, since it does not appear practical to design wrists with more than four revolutes, we 
limit ourselves, in the balance of the paper, to four-revolute manipulators, i.e., we set n = 4, 
unless otherwise stated. 



4 Four- Axis Isotropic Spherical Wrists 

In this section we obtain all possible four-revolute serial spherical wrists with isotropic architec- 
tures. 

The algebraic problem at hand consists in finding the set of vectors { }f that verify the 
isotropy conditions of eq.©. Without loss of generality, we define ei parallel to the x axis of 
the coordinate frame at hand; then, we let e 2 lie in the y-z plane of the same frame, while the 
remaining two vectors are left arbitrary. We thus have 



" 1 " 




c 




X 




u 





. e 2 = 


s 


i e 3 = 


y 


, e 4 = 


V 







_ _ 




z 




w 



ei 



in which c = cosai and s = sinai, as per Fig. [TJ The isotropy condition 
foregoing components and with a 2 — 4/3, yields, then, 



1 + c- + x z + ?r 
s 2 + y 2 + v 2 

2 i 2 
Z + W 

cs + xy + uv 
zy + wv 
xz + uw 

Besides, we have the normality of e 2 and 63: 

2 1 2 
c + s 

2 1 2 1 2 
x + y + z 



4/3 

4/3 

4/3 







1 
1 



(9) 

in terms of the 

(10a) 
(10b) 
(10c) 
(lOd) 
(lOe) 
(lOf) 

(10g) 
(lOh) 



the normality of e 4 being embedded in the foregoing system of equations. Indeed, this is 
obtained upon adding eqs. fTUk -c) and subtracting this sum from the sum of eqs.( |10g| fe h). We 
have now eight quadratic equations for eight unknowns. The Bezout number of this system is 
thus 2 8 = 256, which means that up to 256 solutions are to be expected, including real and 
complex, as well as multiple solutions. Moreover, the BKK bound [7] of the same system turns 
out to be 192. It will be shown presently that this number is too big, the total number of 
solutions being substantially smaller. In order to find the solutions of interest, we eliminate 
successively all the unknowns but u to obtain a monovariate polynomial in this unknown. First, 
we solve for x, y and c from eqs. (|10dl -f). thus obtaining 



z 

vw 



u(wy — vz) 



HZ 



(11a) 
(lib) 

(11c) 
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Upon substituting eqs. pial -c) into eqs. pOal -c) and Q10g|fc b), we obtain 



w 2 u 2 + v 2 w 2 + z A - z 2 = (12a) 
3sV + 3«V + 3dV-4z 2 = (12b) 
3u 2 w 4 v 2 + 6u 2 w 2 v 2 z 2 + 3u 2 v 2 z 4 + 3w 2 u 2 s 2 z 2 + 3u 2 z 4 s 2 - s 2 z 4 = (12c) 

242,o2222i 224,44 24 n 

u w v + 2 u w v z + u v z + s z — s z = (12d) 

3z 2 + 3w 2 -4 = (12e) 

It is noteworthy that the system of equations (fT2"k -e) contains only second and fourth powers 
of all the unknowns, which allows for a recursive solution, as we shall show below. First, from 
cq.JTJfe) we solve for w: 

w = ±^y/l2-9z 2 (13a) 
Upon substitution of eq. (|13ap into eq. p^b ) we obtain s: 



3 z 

Likewise, substitution of eq. (|13bj) into cq. (fT2l i) yields v. 



± 2V3^ M 



v = ±\Jz 2 - 4m 2 (13c) 

Finally, substitution of eq. (|13cp into ea. p^H ) leads to z: 

z = ±VGu (13d) 

Now, substitution of eqs. (|13al -d) into eq. (|12c|) leads to a monovariate polynomial: 

u(3u-l)(3«+l) = (13e) 

The above equation reduces, in fact, to a quadratic equation because u cannot vanish, as we 
shall show presently. Thus, the two possible solutions for eq. (|13ep are 



u = ±l (14) 

We thus have a set of five quadratic expressions for the five unknowns u, z, v, s and w, which 
means that we have found 2 5 = 32 distinct solutions, as displayed in Tabled Therefore, the 
Bezout number of this system overestimates the number of solutions by a factor of eight, while 
the BKK bound by a factor of six. 

The solutions can now be readily computed recursively. Indeed, eqs.([T3a]-c) lead to 

v = ±V2u (15a) 

- = ±^ (15b) 
1 



w = ±- x /6(2-9u 2 ) (15c) 

o 

the remaining unknowns being computed from eqs. (|13d[) . (|lla[) , (|llb[) and (|llc[) . in this order. 
Now we show that none of the unknowns can vanish. We do this by noting that: 



G 



1. If u = 0, then v = by virtue of eq. (|15ap . while w = ±2\/3/3, by virtue of eq. p5cp . 
thereby violating the normality condition on e4. Also note that u = leads to c = 
by virtue of eq. picll . but this value, along with that of s given by eq. (|15bj) . violates the 
normality condition ( |10g[ ); 

2. If v = 0, then u = by virtue of ea. (|15ap . but, according to item[TJ this is impossible; 

3. If w = 0, then, by virtue of ecis. ljllal fc b), x = y = 0; additionally, by virtue of eqs. (|13d[) 
and ([HI) , z = ±\/6/3, thereby violating the normality condition (|10h[k 



Sol'n # 


c 


s 


X 


V 


z 




V 


w 


1 


1/3 


-2V2/3 


-1/3 


-V2/3 


V6/3 


1/3 


y/2/3 


y/G/3 


2 


1/3 


-2V2/3 


-1/3 


-V2/3 


V6/3 


-1/3 


-y/2/3 


-V6/3 


3 


1/3 


-2V2/3 


-1/3 


-V2/3 


-V6/3 


1/3 


V2/3 


-V6/3 


4 


1/3 


-2V2/3 


-1/3 


-V2/3 


-V6/3 


-1/3 


-V2/3 


V6/3 


5 


1/3 


-2V2/3 


1/3 


V2/3 


V6/3 


1/3 


y/2/3 


-V6/3 


6 


1/3 


-2V2/3 


1/3 


V2/3 


V6/3 


-1/3 


-y/2/3 


V6/3 


7 


1/3 


-2^2/3 


1/3 


\/2/3 


-V6/3 


1/3 


V2/3 


V6/3 


8 


1/3 


-2V2/3 


1/3 


\/2/3 


-V6/3 


-1/3 


-V2/3 


-V6/3 


9 


1/3 


2V2/3 


-1/3 


\/2/3 


V6/3 


1/3 


-V2/3 


V6/3 


10 


1/3 


2V2/3 


-1/3 


\/2/3 


V6/3 


-1/3 


V2/3 


-V6/3 


11 


1/3 


2\/2/3 


-1/3 


V2/3 


-V6/3 


1/3 


-V2/3 


-V6/3 


12 


1/3 


2V2/3 


-1/3 


V2/3 


-V6/3 


-1/3 


y/2/3 


V6/3 


13 


1/3 


2V2/3 


1/3 


-V2/3 


V6/3 


1/3 


-y/2/3 


-V6/3 


14 


1/3 


2^2/3 


1/3 


-V2/3 


V6/3 


-1/3 


y/2/3 


V6/3 


15 


1/3 


2^2/3 


1/3 


-V2/3 


-V6/3 


1/3 


-y/2/3 


V6/3 


16 


1/3 


2V2/3 


1/3 


-V2/3 


-V6/3 


-1/3 


y/2/3 


-V6/3 


17 


-1/3 


-2V2/3 


-1/3 


\/2/3 


V6/3 


1/3 


-y/2/3 


V6/3 


18 


-1/3 


-2a/2/3 


-1/3 


V2/3 


\/6/3 


-1/3 


y/2/3 


-V6/3 


19 


-1/3 


-2V2/3 


-1/3 


V2/3 


-V6/3 


-1/3 


y/2/3 


V6/3 


20 


-1/3 


-2V2/3 


-1/3 


V2/3 


-V6/3 


1/3 


-y/2/3 


-V6/3 


21 


-1/3 


-2V2/3 


1/3 


-V2/3 


V6/3 


1/3 


-y/2/3 


-V6/3 


22 


-1/3 


-2^2/3 


1/3 


-V2/3 


V6/3 


-1/3 


y/2/3 


V6/3 


23 


-1/3 


-2V2/3 


1/3 


-V2/3 


-V6/3 


-1/3 


y/2/3 


-V6/3 


24 


-1/3 


-2V2/3 


1/3 


-V2/3 


-V6/3 


1/3 


-y/2/3 


V6/3 


25 


-1/3 


2V2/3 


-1/3 


-V2/3 


-V6/3 


1/3 


y/2/3 


-V6/3 


26 


-1/3 


2^/2/3 


-1/3 


-V2/3 


-V6/3 


-1/3 


-y/2/3 


\/6/3 
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1/3 
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V6/3 


29 
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1/3 


V2/3 


-V6/3 


-1/3 


-y/2/3 


-V6/3 


30 


-1/3 


2^2/3 


1/3 


V2/3 


-V6/3 


1/3 


y/2/3 


V6/3 


31 


-1/3 


2V2/3 


1/3 


y/2/3 


V6/3 


-1/3 


-y/2/3 


V6/3 


32 


-1/3 


2V2/3 


1/3 


y/2/3 


V6/3 


1/3 


y/2/3 


-V6/3 



Table 2: The 32 isotropic sets of four points on the unit sphere 
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4. If x = 0, then, according to eq. ljl3dp . either u = or w — 0, but none of these can vanish, 
according to items Q] and 

5. If y = 0, then, by virtue of eq. ljllbj) either v = or w = 0, but these alternatives are not 
plausible according to items [2] and 

6. If z = 0, then u — 0, according with eq. (|13d[) . but, by virtue of item[ll this is not possible; 

7. If c = 0, then the normality of e2 requires that s = ±1, but this is impossible by virtue of 
eq.(HS9; 

8. From eq. (|15b[) . s cannot vanish, thereby showing that none of the unknowns can vanish. 

Note that the trivial isotropic set of points given in eq. ([2"a[) is solution 18 of the Tabled 
Table [3] records seven isotropic sets of points obtained by means of antipodal exchanges. 





P2 


Pz 


Pi 


P2P3 


P2P4 


P2P4 


P2P3P4 


Solution # 


10 


23 


17 


16 


24 


9 


15 



Table 3: Isotropic sets of points obtained with antipodal exchanges 



With the forgoing eight isotropic sets of points, we obtain additional isotropic sets by reflec- 
tions onto the coordinate planes x-z, x-y, x-z and x-y; we can then verify that these solutions 
are listed in Table [H The corresponding solutions are given in Table [U 



Solution # 


18 


10 


23 


17 


16 


24 


9 


15 


Reflection plane: x-y 


19 


12 


22 


20 


14 


21 


11 


13 


Reflection plane: x-z 


27 


2 


29 


28 


8 


30 


1 


7 


Reflection planes: x-z and x-y 


26 


4 


31 


25 


6 


32 


3 


5 



Table 4: Isotropic sets of points defined by reflections onto the coordinate planes x-z and x-y 

It is noteworthy that Table 4 includes a reflection about the x-y and the x-z planes, which 
amount to a 180° rotation about the x axis, and does not lead to a new wrist. 

Now, for all seven solutions of Table [3] and the trivial isotropic set of points given in eq. (|2"a|) , 
we compute the corresponding Denavit-Hartenberg (DH) parameters yielding isotropic wrists. 
For each isotropic set of points, we place the first joint axis at P\ . It is now apparent that we can 
derive six kinematic chains. Thus, we find cti, ct2 and 013 as the angles made by the neighboring 
position vectors of points Pj. Moreover, we eliminate the set of DH parameters leading to wrists 
that are identical. Hence, a total of eight distinct isotropic wrists are obtained from these sets. 
The Denavit-Hartenberg parameters of the eight distinct wrists are displayed in Table [5] The 
corresponding wrists, at their isotropic postures, being displayed in Figs.[^i-h. 

Note that the entry corresponding to 0:4 in the foregoing table is left with an asterisk because 
this twist angle is not defined for a four-revolute wrist. Its value depends on how the z-axis 
of the task frame is defined. As well, angles B\ and #4 are left unspecified because isotropy is 
independent of these values, i.e., isotropy is preserved upon varying these two angles throughout 
their whole range of values, from to 2-7T. 
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5 Conclusions 



We showed that the algebraic formulation of the problem leading to all four-revolute serial 
spherical wrists with kinetostatic isotropy yields a system of eight quadratic equations in eight 
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Table 5: The Denavit-Hartenberg parameters of the eight spherical wrists at their isotropic 
postures 




(a) (b) (c) (d) 




(e) (f) (g) (h) 



Figure 2: The eight distincts isotropic wrists of Table 5 



9 



unknowns, whose Bezout number is 256, its BKK bound being 192. Nevertheless, this system 
admits only 32 distinct solutions. Furthermore, upon elimination of the solutions leading to 
repeated wrists, we are left with only eight distinct isotropic wrists, whose Denavit-Hartenberg 
parameters were computed and displayed, the corresponding wrists having been displayed at 
their isotropic postures. 
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Appendix 

We show here that a reflection L about a line C that passes 
through the origin is a rotation about C through an angle tt. 
To this end, we resort to Fig. [31 showing C and a point P. To 
simplify matters, we sketch P and C in their plane. 

The projection of P onto C is denoted by P' , the reflection 
sought by P" , the corresponding position vectors being denoted 
by P, p' and p". Apparently, 

p' = ee T p (16) 

Hence, 

p" = p' - (p - ee T p) = (2ee T l)p 

Thus, the reflection sought, L, is given by the matrix coefficient of p in the rightmost side of 
the foregoing equation, i.e., 

L = 2ee T - 1 (17) 

As the reader can readily verify, the above expression yields a proper orthogonal matrix, and 
hence, a rotation. Moreover, the axis of the rotation is e and the angle 7r. 




Figure 3: The reflection of a 
point with respect to a line 
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